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ABSTRACT 

I consider Lagrangians which depend nonlocally in time but in such a 
way that there is no mixing between times differing by more than some finite 
value At. By considering these systems as the limits of ever higher deriva- 
tive theories I obtain a canonical formalism in which the coordinates are 
the dynamical variable from t to t + At. A simple formula for the conjugate 
momenta is derived in the same way. This formalism makes apparent the vir- 
ulent instability of this entire class of nonlocal Lagrangians. As an example, 
the formalism is applied to a nonlocal analog of the harmonic oscillator. 
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1 Introduction 



The traditional goal of fundamental physics is to infer the rules by which 
the "present state" of a system's dynamical variables determines their future 
state. Since Newton's time, most attention has been given to models for 
which the "present state" of a system's dynamical variables means their 
values at some instant in time and possibly also the values of their first 
time derivatives. This restriction corresponds to equations of motion which 
are local in time and contain no more than second time derivatives. It has 
not so far proved useful, in describing the physical universe on the most 
fundamental level, to invoke equations of motion which are either nonlocal 
in time or which even possess more than two time derivatives. 

The deep reason behind this surprising simplification of fundamental the- 
ory seems to be the result obtained by the 19th century physicist Ostrograd- 
ski . He showed that Lagrangians which possess a finite number of higher 
time derivatives, and are not degenerate in the highest one, must give rise to 
Hamiltonians which are linear in essentially half of the canonical variables. 
This is a nonperturbative result. Further, it cannot be altered by quantiza- 
tion since the instability occurs over a large volume of the canonical phase 
space. I will review Ostrogradski's construction in Section 2 of this paper. 
For now it suffices to note that the instability must apply as well to nonlocal 
theories which can be represented as the limits of ever higher derivative ones. 

Much of the interest in nonlocal quantum field theories has been driven 
by the close connection between ultraviolet divergences and local interactions 
M. Of course it does no good to avoid divergences by introducing an infinite 
number of instabilities against which there is not even any barrier to decay! 
It is therefore of interest to know when a nonlocal Lagrangian possesses a 
higher derivative representation and, consequently, the Ostrogradskian insta- 
bility. The higher derivative representation does seem to be valid for cases 
such as string field theory, where the nonlocality enters through entire func- 
tions of the derivative operator and the Lagrangian cannot be made local 
by a field redefinition |3[. On the other hand, the higher derivative repre- 
sentation is certainly not valid for the inverse differential operators which 
result from integrating out a local field variable. It also fails for "maximal 
nonlocality" in which the action is a nonlinear function of local actions [|]]. 
The purpose of this paper is to demonstrate by construction that the higher 
derivative representation is valid for "nonlocality of finite extent" in which 
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the Lagrangian connects no times differing by more than some constant At. 

Although the results of this paper apply as well to field theories I will 
work in the context of a one dimensional, point particle whose position as 
a function of time is q(t). A nonlocal Lagrangian of finite extent At is one 
which definitely depends upon (and mixes) q(t) and q(t+At), and potentially 
depends as well upon q(t') for t < t' < t + At. An example would be the 
following nonlocal generalization of the harmonic oscillator: 

I i 
L[q] (t) = -mq 2 (t + At/2) - -muj 2 q{t)q{t + At) . (1) 

The deterministic way of viewing such theories is that the equations of motion 
give the dynamical variable at the latest time — g(t + At) — as a function of 
earlier times in the range t — At < t' < t + At. In our example the equation 
of motion is 

C dr 5L[< fq\t) r) = ^ + + + ~ At) l = ° ' 

(2) 

and its deterministic interpretation is 

q{t + At) = -q{t-At)-\q\t). (3) 

This paper contains six sections of which this Introduction is the first. 
Section 2 is devoted to a review of Ostrogradski's result for local Lagrangians 
depending upon iV time derivatives. My canonical formalism is presented in 
Section 3 and shown to correctly realize the dynamics of nonlocal Lagrangians 
of finite extent. This formalism is applied in Section 4 to the Lagrangian (j^) 
discussed above. The connection with Ostrogradski's formalism is demon- 
strated in Section 5. My conclusions comprise Section 6. 



2 Ostrogradski's construction 

Consider a Lagrangian L (q,q, . . . , q^j which depends upon the first N 
derivatives of the dynamical variable q(t). I shall assume only that the La- 
grangian is nondegenerate, i.e., that the equation 

dL 

Pn = d^) ' (4) 
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can be inverted to solve for q( N ' as a function of Pn, q and the first N — 1 
derivatives of q. This just means that the action's dependence upon q( N > 
cannot be eliminated by partial integration, so the equation of motion, 



N 



dV dL 



1=0 



contains q^ 2N \ 

Since the equation of motion determines g( 2Ar ) as a function of q and 
its first 2N — 1 derivatives, one can obviously specify the initial values of 
these 2N variables. The canonical phase space must accordingly contain N 
coordinates and N conjugate momenta. In Ostrogradski's construction Jl | 
the 7-th coordinate is just the (I — l)-th derivative of q, 

Qi = ■ (6) 

The momentum canonically conjugate to Qi is, 

A consequence of nondegeneracy is that the derivatives q( N+I ^ can be deter- 
mined from Pn-i, Pjv-j+i> • • • > Pn and the Qj's. In particular q^ involves 
only Pjy and the Q/s, 

q (N) = Q (Q, Pn) ■ (8) 
Ostrogradski's Hamiltonian is, 

N 

H = Y,PiQi~ L i (9) 
1=1 

N-1 

= Y, p iQi+i + p nQ{Q,Pn)-l(q,q(q,p n )) , (10) 

1=1 

and his canonical equations are the ones suggested by the notation, 

6 - dH P - dH (ID 

Ql ~m ' Pl -~dQi- ( } 
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It is straightforward to check that the various canonical evolution equations 
reproduce the equation of motion and the structure of the canonical formal- 
ism: Qj gives the canonical definition (^) for Qi + i, Qn gives the canonical 
definition for P^ in its inverse form (|j); -P7+1 gives the canonical definition 
(0) for Pi] and P\ gives the equation of motion @. So there is no doubt that 
Ostrogradski's Hamiltonian generates time evolution. When the Lagrangian 
is free of explicit time dependence H is also the conserved current associated 
with time translation invariance. 

The instability consequent upon if's linearity in P±, P2, . . . , Pn-i explains 
why higher derivative theories have not been of use in describing physics on 
the fundamental level. Note the generality of the problem. It does not depend 
upon any approximation scheme, nor upon any feature of the Lagrangian 
except nondegeneracy. Further, it must continue to afflict the theory after 
quantization because the instability is not confined to a small region of the 
classical phase space. If a fully nonlocal Lagrangian can be represented as the 
limit of such higher derivative Lagrangians it must inherit their instability. 

The limit of infinite N is facilitated by regarding Ostrogradski's formalism 
as the result of constraining a larger system with an extra pair of canonical 
variables, 

Q N+1 = qW , P N+ i~0. (12) 



The Hamiltonian is, 



TV 

H = J2 P1Q1+1 - L (Q, Qn+i) , (13) 
7=1 

and requiring that P/v+i remains zero imposes the canonical definition of P^ 
as another constraint, 

p ^ = -m^ = - p ^^" a - (14) 

Since the Poisson bracket with Pjv+i gives the second derivative of the La- 
grangian with respect to Qn+i, nondegeneracy implies that the two con- 
straints are second class. The resulting Dirac brackets are, 



{Qi,Qj} d = (Si n+i5j n + 5i N 5 



I NO J N+l) 



d 2 L 



dQ 2 



N+l 



(15) 
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{Qi,Pj} d 
{Pi,Pj} d 



& I J — ^ I N+l 
. 



d 2 L 



N+l 



-1 



d 2 L 



OQjOQ 



N+l 



(16) 
(17) 



Note that there is not even any difference between Dirac brackets and Poisson 
brackets provided one avoids the highest Q — that is, Qn+i- 



3 My construction for finite nonlocality 

I define a nonlocal Lagrangian L[q](t) of finite extent At as one which poten- 
tially depends upon the dynamical variable from time t to time t + At, with 
guaranteed mixing between q(t) and q(t + At). The requirement of mixing 
is the generalization of nondegeneracy and it implies, 

*0. (18) 



5q(t)5q(t + At) 



I shall also require that the Lagrangian contain no derivatives of either q(t) 
or q(t + At). 

I label the canonical variables by a continuum parameter < s < At. 
They are defined as follows: 

Q(s,t) = q(s + t), (19) 

Note that (|20D implies the constraint P(At,t) w 0. Note also that whereas 
the s and t derivatives of Q(s,t) are identical those of P(s,t) are not, 

Since L[q]{t) involves the dynamical variable from q(t) up to q(t + At) we 
see that P(s, t) involves q(s + 1 — At) up to q{t + At). So decreasing s allows 
one to reach back further before time t, all the way to time t — At at s — 0. 

Note that the equation of motion is P(0,t) = 0. This emerges as an 
additional constraint from surface variations of the canonical Hamiltonian, 

f At A 

H(t) = J q drP(r,t)-Q(r,t) - L[Q)(t) . (22) 
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We can find the canonical equations of time evolution from the fact that the 
only nonzero Poisson bracket is, 

{Q(r,t),P(s,t)} = 5(r-s). (23) 

The result for Q(s,t) is straightforward, 

j t Q(s,t) = {Q{s,t),H(t)} , (24) 

rAt A 

= j dr5(r - s)—Q(r,t) , (25) 
= ^QM)- (26) 

A partial integration is necessary for P(s, t) and one must be careful about 
the resulting surface terms, 

j t P(s,t) = {P(s,t),H(t)} , (27) 
'■ M d . SL[Q](t) 



(29) 



For < s < At this simply reproduces ([21]), and hence the canonical defini- 
tion of P(s, t). 

Since the two surface terms cannot be canceled by anything else, they 
must be imposed as constraints, 

P(s,t)&0 s = 0,At. (30) 

Requiring that they be preserved under time evolution implies two additional 
constraints, 

^>+iSH (31) 

Nondegeneracy — and the absence in L[q](t) of derivatives of q(t) and/or 
q(t + At) — guarantees that the four constraints are second class. 
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Note that the H(t) is conserved when L[q}(t) is free of explicit time de- 
pendence, 



dH 
~dt 



(t) 



A * (dP 



dQ, 



d 2 Q 



*<^(M)^(M) + p(M)-^(M) 



d_ 

dt 



L[Q](t)(32) 



A* d 

ds— 

ds 



P(s,t)—Q(s,t) 
ds 



+ 



A * , 8L[Q](t) d , d rr ^ 1/ , 



d 



p(s,t)_g( s ,t) 
as 

. 



Af 



(33) 

(34) 
(35) 

Note also that the Hamiltonian has inherited the Ostrogradskian instability. 
After eliminating the constraints it must be linear in all the P(s,t) except 
possibly 4-P(s,i) at s = and at s — At. 



4 A simple example 



It is useful to see how the general construction given in the previous section 
applies to the Lagrangian ([!]) presented in Section 1. Of course the canonical 
coordinates are always Q(s,t) = q(s + t) for < s < At. To find the 
canonical momenta note that the functional derivative of the Lagrangian is, 



5L[q](t) 
5q(s + t) 



-mq t 



At 



5' 



At\ 



— muj 
2 



q(t + At)S(s) +q(t)5(s- At) 



gives 
At\ 



2 J*' " 



At\ 



At 



mq(s + t)6 y ^ 



Substituting in ( 
P(s,t) = mq\t + 

-Knu 2 q{s + t + At)6(-s) +q{s + t- At)6(At - s) 

Note that P(At,t) = and that 

- i i 

P(0, t) = -m q\t) + -tu 2 q(t + At) + -uj 2 q(t - At) 

2 2 



(36) 



(37) 



(38) 
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indeed vanishes with the equation of motion. 
The canonical Hamiltonian is 



m 

2 



dQ (At 



2 



^mu 2 Q{0,t)Q{At,t) . 



ds V 2 ' 

(39) 

The canonical evolution equations are, 

f(M) = f (.,.). (40) 
dP, x gLP . . dQ [At \ ., / At\ 

-imcj 2 [Q(AM)5(s) + Q(0,t)5(s- At)] . (41) 

It is simple to check that substituting Q(s,t) = q(s + t) and relation Q3"7D for 
P(s, t) indeed verifies these equations. 

The constraints are P(0, t) ~ 0, P( At, i)«0 and the apparently singular 
pair, 

dP 1 

— (0,t) - -muj 2 Q(At,t)5(0) w 0, (42) 
as 2 

dP 1 

— (At,t) - -mcj 2 Q(0,t)o(0) w 0. (43) 

However, the vanishing of P(s, t) at the endpoints means that the endpoint 
derivatives contain delta functions, so the actual constraints are the perfectly 
regular coefficients of 5(0), 

P(0 + ,t) - -muj 2 Q(At,t) w 0, (44) 
-P(At-,t) - -muj 2 Q{0,t) w 0. (45) 

Note that these constraints are implied by ([37]) and, where necessary, the 
vanishing of (|38|). Note also that the contstraints determine both the actual 
endpoint values of P(s, t) and its limit as the endpoints are approached. 

Since the Lagrangian (|l|) has no explicit dependence upon time the Hamil- 
tonian should be conserved. To see that it is, first substitute Q(s, t) = q(s, t) 
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and relation fl3~7|) for P(s,t) to obtain, 



l M dsP(s,t)^-(s,t) 
Jo dq 



-m 



dQ, s -i( ^\ f At/2 , •/ 

' 11 ■■■-'t-\ )-m dsq(s + t)q(s + t) 

2 J Jo 

1 r At 

— muj 2 / dsq(s + t - At)q(s + t) , (46) 

2 Jo 

/■At 

/ rfsg(s + t - At)g(s + t) . (47) 
Jo 



'/ 2 |i + f)+g 2 (t) 



-mio 



Then subtract expression (JTJ) to determine the configuration space Hamilto- 



man, 



11 1 f At 

H(t) = -mq 2 {t) + -muj 2 q(t)q{t + At) - -muj 2 / dsq(s + t - At)q(s + t) . 
2 2 2 Jo 

(48) 

Now use the fact that the integrand depends upon t only through the sum 
s + 1 to express the derivative of the integral as a surface term, 



dH 1 

— (t) = mq(t)q(t) + -muJ 2 {q(t)q(t + At)+q(t)q(t + At)] 



Y s=At 
- -muj 2 q{s + t - At)q(s + t) 



s=0 



mq{t) 



q\t) + ^u 2 q(t - At) + ^u 2 q(t + At) 



(49) 
(50) 



The most straightforward way of demonstrating that the transformation 
to the constrained phase space is invertible is by exhibiting the inverse. Of 
course we always have, 



q(s + t) = Q(s,t) , V0<s<At. 
For —At < s < one recovers q(s + 1) from relation ([37]). 



q(s + t) 



muo 



2 „/ A x 2 d 
—P(s + At,t) 



uj 2 ds 



dQ { 
ds 



s + At,t)9 (~ ~ sj 



(51) 



(52) 



The endpoint case of s = —At is given by the constraint P(0,t) ~ 0, 



q{t-At) = -Q{At,t) 



2 d 2 Q 
uj 2 ds 2 



(53) 



where I am of course defining differentiation in the right-handed sense, 

#(x) = Km /(x + £) ' /(3:) . (54) 

It is amusing to close the section by exhibiting the run-away solutions 
which are one possible consequence of the Ostrogradskian instability. Since 
the configuration space equation of motion, 

q(t) + Kj 2 [q{t + At) + q(t - At)] = , (55) 

is linear and invariant under time translation, the general solution must be 
a superposition of terms having the form e lkt . The allowed frequencies are 
complex numbers k which obey, 

k 2 = u 2 cos(kAt) . (56) 

The equation is transcendental but graphing both sides shows a single pair 
of ± real solutions. To find the remaining solutions make the substitution, 

k = a + i(3 , (57) 

and take the real and imaginary parts of the equation, 

a 2 - (3 2 = oj 2 cos(aAt) cosh(/?At) , (58) 
2a(3 = -oj 2 sin(aAt) sinh^At) . (59) 

Graphical analysis indicates a conjugate pair of solutions for a At in each 2tt 
interval of the real line. For large integer N these solutions have the form, 

21n(iV) , . 

a At w 2txN , 60 

ttN 

, f87i 2 N 2 \ f\n(N)\ 2 

^ At ^ ln (^AFj + (^rj • (61) 

So this system has the infinite number of solutions predicted by the Ostro- 
gradskian analysis, and all but two them grow or fall exponentially. 
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5 Ostrogradskian derivation 

My representation is related to the infinite N limit of Ostrogradski's through 
the Maclaurin series, 

oo I 

Q(s,t) = J2-»Qi+i(t)- (62) 
7=0 1 - 

Note that differentiation with respect to the Ostrogradskian coordinates is 
realized by the functional chain rule, 



dQ^f) ~ ' "' S 



dQ(s,t) 
dQi(t) 



ds— — — r^— — , (63) 



5Q(s,t) Jo (I-l)\6Q(s,t) 



where the functional derivative is defined by 

5Q(r,t) 



SQ(s,t) 



5{r - s) , (64) 



and the ordinary rules of calculus. From (|i~9|) one obtains a useful formula 
for the higher derivative representation, 

am*) _ t*> ds °' «MM (65) 



dqW(t) Jo I\5q(s + t) 

The conjugate momentum P(s, t) should depend linearly on the Ostro- 
gradskian momenta, 

oo 

P(s,t) = y £p I (8)Pi +1 {t). (66) 

7=0 

The combination coefficients pi(s) can be determined by enforcing the canon- 
ical Poisson bracket (EBl), 



oo J oo 

7=0 1 ■ J=0 

oo r I 

= Enf^)- (68) 
7=0 1 ■ 



By acting (d/dr) J and then taking r — > one finds, 
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To obtain my formula ( p0|) for the conjugate momenta note first that, for 
infinite N, the Ostrogradskian momenta are, 

oo I J— I 

Plit) 



00 / d 
hi \dt 



dL[q)(t) 



Now substitute this and fl69| ) into (|6 



dq( J \t) ' 
J -< -^ dr r J SL[q](t) 



J\ 5q(r + t) 



1=0 



I) 



d 



J-I-l 



s [-it 

j=i+i \ UL 



At r J 



dr 



r SL[q](t) 
J\ 5q(r + t) 



Simplification is achieved by exploiting the identity, 



T\ J,, 



dr 



,(r-r')V) 



7!(J-/-1)! ' 

to recognize the two sums as Taylor expansions of the shift operator, 

i-At 

lo 



P(s,t) 



dr / dr' ^ 

j=o 







00 

x E 

j =1+1 



(r — r') 1 


( da) 


I\ 




(r>y~i- 


u- 


(J-I- 





S(s) 



dt 



j-7-1 



mm 

5q(r + t) 



dr dr'5{s -r + r') - / n 

Jo oq(r +t — r) 



A * , 8L[q](s + t-r) 

dr — r . 

5q{s + t) 

The Hamiltonian follows similarly, 

Hit) = Y, p MQim-L[Q]{t) 
1=1 

00 r At J-1 / 



E 

1=1 

i-At 



ds 



d 



d 



P(s,t)l- \ Q(r,t) 



r=0 



r m d 

J q dsP(s,t)—Q(s,t) - L[Q}(t) . 



(70) 
(71) 

(72) 
(73) 



(74) 

(75) 
(76) 

(77) 

-L[Q](t), (78) 
(79) 



Its instability is manifest from the fact that it has been derived from Ostro- 
gradski's result in the limit that the number of derivatives becomes infinite. 
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6 Discussion 



I have shown that Lagrangians with nonlocality of finite extent At can be 
treated as the limits of higher derivative Lagrangians. I have also given a 
canonical formalism that is somewhat more natural in which the canonical 
variables are labelled by a continuum parameter s, for < s < At. The 
canonical coordinates are just the dynamical variables at times t + s. A 
quantum mechanical state in such a system would be a functional of these 
coordinates. The conjugate momenta Q20D are given by a simple integral of a 
functional derivative of the Lagrangian. With the canonical coordinates, the 
momenta allow one to reconstruct the dynamical variables at times t — s. 

There is no physical motivation for this exercise because all such models 
are virulently unstable. Indeed, the only point of the formalism is to remove 
any doubt about a possible phenomenological role for these Lagrangians. 
They have inherited the full Ostrogradskian instability: essentially half of 
the directions in the classical phase space access arbitarily negative energies. 
There is not even any barrier to decay. This is a nonperturbative result 
and, because it arises from a large region of phase space, it must survive 
quantization. 

Negative results of such power and generality are to mathematically in- 
clined physicists like a red flag is supposed to be to a bull. Nothing I can 
honestly add is likely to much discourage further attempts to carve out a 
physical niche for nonlocal Lagrangians but I do recommend that these ef- 
forts be preceded by sober reflection upon the following fact: in the long 
struggle of our species to understand the universe it has never once proven 
useful to invoke a theory which is nonlocal on the most fundamental level. 
Yet the subset of local Lagrangians containing no more than first derivatives 
is a minuscule fragment of the set of all functionals of the dynamical variable. 
The Ostrogradskian instability offers a simple and compelling explanation for 
the complete dominance of this tiny subset over its much larger whole. The 
only alternative would seem to be coincidence on a scale that makes even the 
worst fine tuning problem seem inconsequential. 

Note Added: Shortly after posting this paper I learned of important work 
by Llosa and Vives || on the problem of canonically formulating a general 
nonlocal Lagrangian. My work can be viewed as a specialization of their 
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technique to the case of nonlocality of finite extent where the Euler-Lagrange 
equations are deterministic, where an explicit Poisson bracket structure can 
be determined and where the formalism can be derived from the infinite N 
limit of Ostrogradski's construction. (None of these features can be present 
in the general case.) Note should also be taken of the recent work of Gomis, 
Kamimura and Llosa on canonically formulating spacetime noncommutative 
theories ||. 
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